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We present the results from an optical spectroscopic analysis of the mas- 
sive stars HD 37366 and HD 54662. We find that HD 37366 is a double-lined 
spectroscopic binary with a period of 31.8187 ± 0.0004 days, and HD 54662 is 
also a double lined binary with a much longer period of 557.8 ± 0.3 days. The 
primary of HD 37366 is classified as 09.5 V, and it contributes approximately 
two-thirds of the optical flux. The less luminous secondary is a broad-lined, early 
B-type main-sequence star. Tomographic reconstruction of the individual spec- 
tra of HD 37366 reveals absorption lines present in each component, enabling 
us to constrain the nature of the secondary and physical characteristics of both 
stars. Tomographic reconstruction was not possible for HD 54662; however, we 
do present mean spectra from our observations that show that the secondary 
component is approximately half as bright as the primary. The observed spectral 
energy distributions (SEDs) were fit with model SEDs and galactic reddening 
curves to determine the angular sizes of the stars. By assuming radii appropriate 
for their classifications, we determine distance ranges of 1.4 - 1.9 and 1.2 - 1.5 
kpc for HD 37366 and HD 54662, respectively. 

Subject headings: binaries: spectroscopic — stars: early-type — stars: individual 
(HD 37366, HD 54662) 



1. Introduction 

There remains considerable uncertainty about the masses of the most massive stars 
because of the relati vely small number of known binary systems for which accurate masses 



can be determined (jGiesI 120031 ) . Spectroscopic measurements alone yield mass functions 
dependent on the unknown orbital inclination, and the determination of inclination requires 
either the good fortune of finding eclipsing binaries or the angular resolution of the orbit on 
the sky. The angular semimajor axis of a binary (in units of milliarcsec = mas) is given by 

/ s noo Mig d)) 2/3 (M total /(3OM )) 1 / 3 

a(mas) = °- 28 m^j) (1) 

where P is the orbital period, M to tai is the combined mass of the stars, and d is the dis- 
tance. The denominators of each unit give typical values for these parameters among OB 
binaries, and the leading coefficient of 0.28 mas indicates that most massive systems are 
probably too closely separated for direct resolution with optical long baseline interferome- 
ters, where the limits are currently above 1 mas. The key objective here is to find double- 
lined spectroscopic binaries with long orbital periods. Such binaries are difficult to detect 
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because their orbital semiamplitudes are small, the component lines are often blended, and 
a long term observational program is required to obtain ad equate phase cov erage. The 



25 



best candid ates for direct resolution are 15 Mon (P 

(P = 442 d; iGarmanv fc Massevlll98ll : be Becker et al-lbood ). and HD 193322 (P 



McKibben et al. 



Gies et al.lll997f l. HD 15558 

311 d; 



1998 



Here we report on new orbits for two such long period massive bi naries, HD 37 366 



and HD 54662. The star HD 37366 (BD+30°968; HIP 26611; 09.5 V, IWalborn 



1973T) is 



a me mber of the Aur OBI association at a distance of approximately 1.3 kpc ([Humphreys 
19781 ). This association has many bright early-type giants and supergiants, but HD 37366 
has the earliest spectral type a mong the member stars that still reside on the main se- 
quence. The Hipparcos mission (lESAlll997l ) detected a visual companion to H D 37366 with 



AH p = 3.5, a separation of 0V58, and a period of approximately 1300 years (IMason et al. 



998 ). The brig 



iter component of these two stars (H p = 7.7) is a radial velocity variable 



Petrie fc Pearce!ll96ll:lYoung|ll942h. and it is known to show asymmetry in its spectral lines 



( iGrigsby. Morrison, fc Andersonlll992l ). Observations with the International Ultraviole t Ex- 



plorer (IV E) confirm that the spectrum displays double lines (jStickland Lloydll200ll ) 



The second target HD 54662 (BD-10°1892; HIP 34536; LS 197; 06.5 V. IWalbornlll972[ ) 
is also the brightest and earl iest member of its resident association, CMa OBI, at a distance 
of 1 . 3 kp c ( Humphrevslll978n. Radial velocity measurements for HD 54662 extend back many 



decades (IPlaskett 



1924 



Conti. Leep. fc Lorrdll977l; iGarmanv. Conti. &: Masseyill980l ). and 



these display only modest variability. However, iFullertonl (Il990i ) noted the presence of blue 
extensions to the spectral lines that probably indicate the presence of a companion in a 
l ong period orbit. The s catter in the IUE velocities also indicates that the star is a binary 
rtstickland fc LloydlboOlh . 



Here we present an analysis of the radial velocities and spectra of both stars from 
spectroscopic observations that we have obtained over the past few years (§2). We discuss 
each system's orbital velocity solution (§3) and the spectral and physical properties of each 
component star in these binaries (§4,5). We conclude with a consideration of the prospects 
for the angular resolution of the orbits using optical long baseline interferometry (§6). 



2. Observations 

We observed HD 37366 and HD 54662 with the Kitt Peak National Observatory (KPNO) 
0.9 m Coude Feed telescope during two separate observing runs in 2000 October and 2000 
December. The spectra were made using the long collimator, grating B (in second order with 
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order sorting filter OG 550), camera 5, and the F3KB CCD, a Ford Aerospace 3072 x 1024 
device with 15/im square pixels. The setup yielded a resolving power of R = A/5A=9500, 
with a spectral coverage of 6440—7105 A. Exposure times were usually 10 minutes or less, 
and we generally obtained two spectra (taken a few hours apart) each night. For HD 37366, 
we made two more red spectral observations in 2004 October using a similar arrangement 
but with a different detector, the T2KB CCD (2048 x 2048 24/im square pixels). In 2006 
October, both HD 37366 and HD 54662 were observed in the red region again using this 
same instrumental set-up. We also observed the rapidly rotating A-type star, ( Aql, which 
we used for removal of atmospheric water vapor and O2 bands. Each set of observations was 
accompanied by numerous bias, flat field, and Th-Ar comparison lamp calibration frames. 

We also obtained a small set of blue spectra of these targets. For HD 37366, the first 
group of four spectra were made in 2005 October with the KPNO 2.1 m telescope and 
GoldCam spectrometer. We used the #47 grating in second order, recording the spectral 
region from 4050— 4950A with a resolving power of R = X/8X ~ 3000. Then in 2005 
November and 2006 October we obtained higher resolution observations in the blue with the 
KPNO Coude Feed 0.9 m telescope. HD 37366 was observed on both occasions, whereas 
HD 54662 was only included during the 2006 observing run. We used grating A in second 
order with order sorting filter 4-96, camera 5, and the T2KB CCD. This setup gave us a 
resolving power of R = X/5X ~ 12100 and a wavelength coverage of 4240-4585A. 

The spectra were extracted and calibrated using standard routines in IRAF@. All the 
spectra were rectified to a unit continuum by fitting line-free regions. The removal of atmo- 
spheric lines from the red spectra was done by creating a library of ( Aql spectra from each 
run, removing the broad stellar features from these, and then dividing each target spectrum 
by the modified atmospheric spectrum that most closely matched the target spectrum in a 
selected region dominated by atmospheric absorptions. The spectra from each run were then 
transformed to a common heliocentric wavelength grid. 



6 IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Associa- 
tion of Universities for Research in Astronomy, Inc., under cooperative agreement with the National Science 
Foundation. 



- 5 - 



3. Radial Velocities and Orbital Elements 
3.1. HD 37366 

We measured radial velocities of the high resolut ion red spectra collected in 2000 and 



2006 using a template- fitting scheme (IGies et al.ll2002l ) for the He I A6678 line. We decided 



not to measure the other strong lines in this region because the binary components are 
badly blended in the Ha profile and the He I A7065 line was marred by residual features left 
behind by the telluric cleaning procedure. This radial velocity measurement scheme assigns 
template spectra that are approximate matches for the primary (hotter and more massive 
star) and secondary spectra, and then makes a non-linear least-squares fit of the shifts for 
each component that best match the observed line profile. We need to make assumptions at 
the outset about the temperature, gravity, projected rotational velocity, and flux contribution 
of each star, but these can be checked after completion of the velocity analysis by studying 
the properties of tomographically reconstructed spectra of the components (§4). 

The matching template spectra for the primary and sec ondary componen t s wer e con- 



structed from the grid of O-type star model spectra from lLanz fc Hubenyl (120031 ) that 



are based upon the line blanketed, non-LTE, plane-p arallel, hydrostatic atmosphere code 



TL USTY and the radiative transfer code SYNSPEC ( Hubenyl 11988 : iHubenv & Landll995 



Hubeny. Heap, fc Lanall998l ). We selected the spectrum taken on HJD 2,451,901.92, which 
shows well separated, individual components of each star, as a reference to determine the 
approximate spectral parameters for both stars. 

The template fitting procedure also requires preliminary estimates of the primary and 
secondary stars' radial velocities. We estimated these for each spectrum with well separated 
lines using the IRAF splot routine and deblend option to fit two Gaussians to each composite 
profile. We also measured relative radial velocity shifts of the strong interstellar lines in all 
the spectra referenced to the first spectrum in the stack. We then used these relative shifts in 
the interstellar lines (which should remain motionless) to make additional small corrections 
for the wavelength calibrations (all these corrections were < 2 km s _1 ). 

The final radial velocities from this template fitting procedure (the majority of the ob- 
servations) are listed in Table 1 along with the heliocentric Julian date of mid-observation, 
the corresponding orbital phase, and the residual from the orbital fit (observed minus cal- 
culated) for both the primary and the secondary. The typical errors in these velocities are 
also listed in Table 1. We measure only one line for this data set, so we list the charac- 
teristic errors (not individual errors), which are based upon the scatter in closely spaced 
pairs of observations. These errors are 1.3 and 2.2 km s _1 for the primary and secondary, 
respectively. 
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This template fitting routine was also used in determining radial velocities for the high 
resolution blue spectra (collected in 2005 November and 2006 October), using the four lines 
O II A4349, He I AA4387,4471, and Mg II A4481. We followed the same procedure in 
obtaining the spectral templates and the preliminary radial velocity estimates as described 
above. Since no strong interstellar features are apparent in this region, no additional radial 
velocity correction was applied. The line-to-line la errors in these Vr measurements are 
< 1 km s _1 for the primary, and 4 — 9 km s _1 for the secondary (Table 1). These final Vr 
measurements are also presented in Table 1. 

The four observations made in the blue during 2005 October had a much lower reso- 
lution, thus making it difficult to apply this method of template fitting. To avoid possible 
errors from unseen line blending in the two components, we chose to measure only radial 
velocities of the He II AA4541, 4686 lines present since these lines are found only in the spec- 
trum of the much hotter, primary star (§4). We used a parabolic fitting routine to determine 
the mean velocities of these lines (Table 1). The line-to-line lex errors associated with these 
measurements are < 6 km s _1 (Table 1). 

The two red observations made in 2004 showed no indication of double-lined profiles. 
In this case, we measured velocities only for the primary star by parabolic fitting of the 
line cores of He I AA6678, 7065, in order to minimize the influence of the secondary on the 
line profile. The line-to-line la error associated with these fits are < 2 km s _1 (exclusive of 
blending errors). These velocities are also presented in Table 1. 

The final two spectra of HD 37366 were collected and downloaded from the archive of 
the International Ultraviolet Explorer (IUE) satellit^l We measured radial velocities for 
these tw o high dispersion, s hort wavelength, prime camera spectra using a cross-correlation 



method (jPenny et al.lll999l ) with the spectrum of HD 34078 as the reference template. The 
spectrum was double-lined in the first spectrum, SWP 30165. The errors are approximately 
5 km s -1 for the primary and 10 km s _1 for the secondary. These final velocities are also 
presented in Table 1. 

The radial velocities from all the data sets (six total) span 20 years with 45 radial 
velocity measurements for the primary and 38 radial velocity measurements for the secondary 
(Table 1). We first constructed a power spectrum using all the primary star's radial velocity 
measurements, being more reliable and plentiful, to identify possible orbital periods for 
the binary. We used the discret e Fourier transform and CLEAN deconvolution algorithm 



( jRoberts. Lehar. fc Dreherll 19871 ). which shows that the strongest signal occurs near P = 31.7 



days. We then used this estimate as a starting value for the period in fits of the orbital 



7 http:/ /archive. stsci.edu/iue/ 
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elements. 



We determined the orbi tal elements of the binary u sing the non-linear, least-squares, 



orbital fitting program from iMorbey &: Brosterhud (I1974I ). We began with a fit of the pri- 



mary's velocities that is given in column 2 of Table 2 done with equal weighting except for 
the low resolution blue spectra, which have a weight set to zero. This solution has a period 
of 31.8187 ±0.0004 days. The independent orbital solution for the secondary has a period of 
31.822 ±0.002 days, given in column 3 of Table 2. Since the independent solutions agree well 
with each other, we derive a joint solution by fixing the weighted means of the shared orbital 
parameters (P, T, e, uj) found in the independent solutions for the binary in order to make 
fits of the systemic velocity, 71,2, and the semiamplitude, Ki^, for each component (column 4 
of Table 2). In the case of massive binaries, the systemic velocities of the components may 
not agree exactly because of differences in their expanding atmospheres and/or in our case, 
differences in the shapes of the template spectra for the He I A6678, He II A6683 blend. 
The radial velocity curves for the joint solution are plotted together with the observations 
in Figure 1. We also made similar fits weighting each each point by the normalized, inverse 
square of its associated error, these results matched within errors of those from the equal 
weighting fits given in Table 2. 



3.2. HD 54662 

We obtained relative radial velocities for the primary (hotter, more massive) star in 
HD 54662 by cross-correlation with a single spectrum of the star that had good S/N prop- 
erties. These relative velocities were transformed to an absolute velocity scale by adding the 
mean velocity measured through parabolic fits to the cores of the absorption lines in this 
reference spectrum. All the strong lines were included in the cross-correlation measurements, 
namely, Ha, He I A6678 + He II A6683, and He I A7065. We excluded He II AA6527, 6890 
because their measurements deviated from the set listed above, as well as compared to 
each other. We suspect that the residual telluric lines in the spectra, that are very promi- 
nent in these regions, are the cause of this disagreement. The velocities for the primary 
star are presented in Table 3, along with the average velocity and a ( line-to-line) from the 



C IV AA5801, 5812 and He I A5876 lines presented bv iFullertonl (1199(1 ) . We were unable to 
measure velocities for the secondary star in individual spectra due to severe line blending 
with profiles of the primary star (see §5). 

Published velocities for HD 54662 (§1) do not show significant variations. Tab le 3 shows 



that o ur measurements change only slightly over our observation period. However, iFullerton 



(jl990l ) found convincing evidence that this system is a double-lined binary with either a long 
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period or high eccentricity, since he observed a blue-shifted secondary component (suspected 
07 spectral type) in the profiles of C IV AA5801, 5812 and He I A5876. 

Here we present a preliminary orbital solution for the primary compone nt that was de- 



termined using our measurements combined with published m easurements (jPlaskettl 11924 



Garmany et al.l Il980l ; iFullerton 



1990 



Stickland fc Llovdl l200lh for a total of 67 radial ve- 
locities spanning 85 years. IStickland &: Llovdl ((20011) proposed a possible period of ~ 92 d, 



however, their orbit was determined excluding selected data points. We re-investigated the 
possible period by power spectrum analysis of all the available data. We examined all the 
peaks in the CLEANed spectrum using the non-linear, least-squares, orbital fitting routine, 
and among the periods limited by the timescales sampled in our two long runs, we find 
that the best so lution occurs at a period of ~558 days. This confirms the suggestion from 
Fullertonl (119901 ) that HD 54662 is in fact a long period binary. Table 4 lists the preliminary 
orbital elements for HD 54662 assuming equal weighting for all velocities, and this solution 
is plotted in Figure 2. We show below (§5) that these results are affected by line blend- 
ing, and the derived semiamplitude, for example, is a lower limit to the actual value. It is 
also possible that the results collected in the literature have systematic differences related 
to the specific lines and measurement techniques used. These systematic offsets are likely 
much smaller than the system semi-amplitude, and since this system has such a long orbital 
period, we include all available measurements for this preliminary orbital solution. 



4. Tomographic Spectral Reconstruction and Stellar Parameters for HD 37366 



We used a Doppler tomography algorithm (IBagnuolo et al.lll994j ) to separate the pri- 
mary and secondary spectra of HD 37366. We applied tomographic reconstruction to the 
red spectra collected in 2000 (30 total) and to the high dispersion blue spectra collected 
in 2005 and 2006 (five total). Figure 3 shows the reconstructed red spectra for the pri- 
mary (top) and the secondary (bottom). The region affected by the atmospheric band from 
~ 6850 — 7000 A was set to unity. The secondary spectrum shows the weak lines of O II 
AA6641,6721, and C II AA6578, 6582. These lines are absent in the primary spectrum that 
shows instead features such as He II A6683 that are found in O-type spectra. To determine 
a monochromatic flux ratio, F 2 /Fi, we used the equivalent width of He I A6678, sinc e it does 



not change significantly with spectral type for late O- to early B- stars ( 1Contilll974l ). These 
equivalent widths in the primary and secondary reconstructed spectra are equal for a flux 
ratio of F 2 /F 1 = 0.35 ± 0.05. 

We fit these reconstructed spectra with the TLUSTY/SYNSPEC model synthetic spec- 
tra (see §3) to estimate the projected rotational velocity Vsmi, effective temperature T eff , 
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and gravity log g. These values are listed for both components of HD 37366 in Table 5 (where 
subscript 1 identifies the primary and 2 the secondary). For stars like these, the disappear- 
ance of the C II and O II lines and the emergence of the He II and Si IV lines with increasing 
temperature provide a useful temperature estimate, while the width of the Ha wings is sen- 
sitive to the adopted gravity. Vsmi was measured using a rotational broadening function 
applied to the model spectra to fit the two He I absorption lines. The red spectra were first 
used in the determination of these parameters, and the results were later checked with the 
reconstructed blue spectra, which include H7 as well other lines from heavier elements. The 
small Vsini es t imate we derive for the primary agr ees with the IUE measurements from 
Howarth et al.l (119971 ) and IStickland &: Lloydl (1200 ll ) and is much smaller than the value 
for the broader-lined secondary. The primary' s temp erature is somewhat larger than the 
T eff = 29.0 ±1.8 kK estimate by iGrigsby et al.l (J1992J) but the gravities agree exactly. Our 
results for T cS and logg using the TLUSTY c ode are expected to be more reliable than the 
previous models used in IGrigsby et al.l (119921 ). which used the PAM code ( ]Andersonlll985l ) 
that only includes nine elements and many fewer metal lines than does TLUSTY. 

The reconstructions from the five high-resolution blue spectra are presented in Figure 4 
along with identifications of absorption lines. The secondary spectrum (bottom) has lower 
S/N, but even with only five spectra the tomography algorithm was able extract its spectrum. 
It is again apparent that the lines of the secondary are much broader than those of the 
primary. Notice also the absence of the He II A4541 line in the secondary's spectrum, 
reinforcing our conclusion that the secondary is the cooler of the two stars. Based upon the 
secondary's cooler temperature and high surface gravity, we estimate that it is a BO— 1 V 
star. Note that the magnitude difference we derive is larger than expected f or main sequence 
stars separated by only a subtype or so (IMartins. Schaerer. fc Hillierl 120051 ) . so it is possible 
that the primary is a somewhat evolved, more luminous star and/or the companion is a very 
young star close to the ZAMS. It is interesting to note that the high Mi sin 3 i and M 2 sin 3 i 
values from the orbital solution suggest that the inclination is large, % = 60 — 90°. However, 
Hipparcos photometry plotted with the period from our spectroscopic orbital solution shows 
no evidence of eclipses. 



5. Stellar Parameters for HD 54662 from Composite Profile Fits 



Radial velocities measured for HD 54662 were use d to create mea n spectra for our 
observations made in 2000 (Fig. 5) and for those made by iFullertonl (119901 ) in 1986. Figure 6 
shows an expanded view of the regions surrounding the He I profiles for two epochs of 
observation. We see that the secondary component appeared blue-shifted during the 1986 
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run (left panel) and red-shifted in recent spectra (right panel). 

We made preli minary two com ponent fits of the blended He I lines (A5876 for the 
spectra obtained by iFullertonl Il99pl and A7065 for this wor k) using TLUSTY/S YNSPEC 
models. We used the temperature and gravity calibrations of iMartins et al.l (120051 ) to select 
parameters for the composite model profiles to fit our observations. Our model s pectra 
for the primary star are based upon an assumed type of 06.5 V ( Walbornl Il972l ). We 
constructed model spectra for the secondary for spectral subtypes of 07 V - 09.5 V. Next, 
we compared our o bserved mean line p rofiles to these models applying the appropriate flux 
ratio (from AMy in lMartins et al. 1 120051 ) for each spectral component in the shifted, combined 
line profiles. In each trial for a given secondary spectral type, the only variables were the 
component radial velocities and the secondary's p rojected rotational velocity V2 sin i (we 
assumed V\ sin i = 70 km s 1 : ICoriti fc Ebbetall977l ) . Our best match for the secondary was 
made with an 09 V subtype and V^sini = 110 ± 10 km s _1 , which yields a flux ratio of 
F2/F1 = 0.51. Our fits of He I A7065 required us to make small and equal adjustments to the 
model line depths. The resulting fits are shown in Figure 6. We caution that an uncertainty 
in V^sini has a large effect upon the best fit line shifts and flux ratio results. 

The wavelength shifts made to fit these composite line profiles provide us with average 
velocities for the primary and secondary components for each observing run. Assuming that 
the true anomaly v and the longitude of periastron u are known from the preliminary orbital 
fit (Table 4), we may estimate the systematic velocity 7 and semiamplitude K by making a 
least-squares, linear fit of these three velocities using 



V r = 71,2 ± Ki >2 [cos(i> + to) + e cos u] 



(2) 



This solution gives semiamplitudes of K\ = 29 ± 4 km s _1 and K 2 = 75 ± 7 km s _1 and 
systemic velocities of 71 = 45 ± 3 km s _1 and 72 = 40 ± 6 km s _1 . This estimate of the 
secondary radial velocity curve also allows us to compute the component minimum masses 
of the system, Misin 3 ^ « 41.5 ± 7.6 M© and M 2 sin 3 i « 16.0 ± 3.4 M . The radial 
velocity curves for these solutions for the primary (dashed line) and secondary (dot-dashed 
line) are also plotted in Figure 2, along with the time-averaged radial velocities from the 
two-component fits. This analysis of the line blending problem clearly illustrates how the 
presence of the blended secondary spectrum skews the velocity measurements for the primary 
(Table 3) towards the system's center of mass, resulting in a semiamplitude (Table 4) that 
is approximately a factor of two smaller than the actual value. 
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Discussion 



One of the motivations for this study was to find long period binaries that may be 
resolved by optical long baseline interferometry. The CHAR A Array, for example, can re - 
solve binaries with angular separations as small as 1 mas (jten Brummelaar et al.l 120051 ). 
To determine the angular separation of the binaries components, we re-estimated their dis- 
tances by fitting their observed spectral energy distribution (SED) with a model SED to find 
the angular stellar diameters that we then compared with stellar radii estimates for their 
spectral classifications. For each binary, the model temperatures, gravities, and flux ratios 
were applied to create a combined mo del flux dis t ributio n over a range of 1200 — 30000 
A. The galactic extinction curve from iFitzpatrickl ( 19991 ) was then applied to the model 



SED to fit the observ ed photometry for each target. Th e observed SED includes u 



2006 



Cutri et al 



Hauck fc MermiUiodl 119981 ). and 2MASS JHK infrared magnitudes (ISkrutskie et al. 



jravio- 



let flu xes (IUE\ TD-1. iThompson et al.lll978l) and UBV flNeckel. Klare. fc Sarcander 
uvby 



1980) 



2003), all of which were transformed into calibrated flux measurements 



( jColina. Bohlin. fc Castellilll996l : lQravill998l : ICohen. Wheaton. fc Megeathl 120031 ). The best 
fit parameters for reddening E(B — V), ratio of total-to-selective extinction R, and the limb 
darkened angular diameter for the primary Old (from the flux normalization) are listed in 
Table 6. Figures 7 and 8 show the SED plots of these best fits for HD 37366 and HD 54662, 
respectively. 

We then compared the expected theoretical radii for the primar y stars, based upon their 
spectral classifications as main sequence stars ([Martins et al.l 120051 ) . with the angular sizes 
to obtain distances of 1.38 and 1.23 kpc for HD 37366 and HD 54662, respectively (Ta- 
ble 6). These estim ates are consistent with the accepted distances to their home associations 
( [Humphreys! 119781 ) . The distance ranges given in Table 6 reflect the change in stellar radius 
between ZAMS luminosity class V and giant luminosity class III. The mid-point increase in 
size is AR w 3R & for our stars, and we adopted this difference to estimate the associated 
range in distance. Errors will also result from the spread in radius for each spectral sub-type 
bin, but this is quite small (±0.4i? Q ) compared to the luminosity range. The error in 6 LD 
from the SED fit contributes only « 1% in the distance error budget. 

The binary semimajor axis a was found usin g Kepler's Thi r d Law , the derived or- 
bital period, and the stellar mass calibrations from IMartins et al.l (120051 ) (for O stars) and 
Harmanec (|l988h (for B stars). The results for the maximum angular separation p max for 
the projected elliptical orbit are also presented in Table 6 where we give the range in p max 
associated with the range in distance. These separations are too small for speckle resolution 
(p > 0'/035 for Am < 3.0), but they are close to or above the limits of long baseline inter- 
ferometry. The HD 54662 binary system in particular may prove to be important target for 
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mass determination by interferometry. 
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Table 1. HD 37366 Radial Velocity Measurements 



HJD 

( — 2,400,000) 


Telescope/ 
Band 


Orbital 

Phase 


Vi 

(km s _1 ) 


cri 

(km s _1 ) 


(O - C)i 
(km s _1 ) 


Vi 

(km s~ 1 ) 


CT2 

(km s~ x ) 


(O - C) 
(km s _1 


46821.612 


IUE/UV 


0.303 


76.9 


5.0 


0.7 


—57.3 


10.0 


4.3 


46866.116 


IUE/UV 


0.702 


13.2 


5.0 


—0.5 








51817.934 


CF /Red 


0.327 


75.9 


1.3 


— 1.3 


—64.5 


2.2 


— 1.5 


51818.938 


CF /Red 


0.358 


76.2 


1.3 


— 1.0 


—67.1 


2.2 


—4.1 


51819.929 


CF /Red 


0.389 


75.5 


1.3 


-0.6 


—53.2 


2.2 


8.3 


51820.922 


CF /Red 


0.420 


74.7 


1.3 


0.8 


—63.8 


2.2 


—5.2 


51821.918 


CF /Red 


0.452 


71.2 


1.3 


0.4 


—51.9 


2.2 


2.5 


51822.918 


CF /Red 


0.483 


64.1 


1.3 


-2.7 


—49.4 


2.2 


—0.3 


51823.853 


CF /Red 


0.513 


62.2 


1.3 


—0.1 


—38.1 


2.2 


4.9 


51823.980 


CF /Red 


0.517 


60.1 


1.3 


— 1.5 


—41.7 


2.2 


0.4 


51824.881 


CF /Red 


0.545 


55.8 


1.3 


-0.7 


—38.8 


2.2 


—3.6 


51824.997 


CF /Red 


0.549 


57.1 


1.3 


1.3 


—36.9 


2.2 


—2.6 


51830.904 


CF /Red 


0.734 


0.2 


1.3 


— 1.1 


45.9 


2.2 


7.4 


51889.881 


CF /Red 


0.588 


47.9 


1.3 


0.5 


—26.0 


2.2 


—3.0 


51890.819 


CF /Red 


0.617 


42.0 


1.3 


2.0 


—20.9 


2.2 


—7.6 


51892.787 


CF /Red 


0.679 


20.9 


1.3 


—0.8 


17.9 


2.2 


6.7 


51893.855 


CF /Red 


0.713 


9.6 


1.3 


-0.2 


26.5 


2.2 


-0.7 


51894.780 


CF /Red 


0.742 


— 1.6 


1.3 


0.3 


53.6 


2.2 


10.9 


51894.856 


CF /Red 


0.744 


—3.3 


1.3 


—0.5 


47.2 


2.2 


3.2 


51895.874 


CF /Red 


0.776 


-19.6 


1.3 


-2.4 


58.7 


2.2 


-4.5 


51896.802 


CF /Red 


0.805 


—32.7 


1.3 


— 1.0 


78.4 


2.2 


—4.2 


51896.914 


CF/Red 


0.809 


—34.6 


1.3 


— 1.0 


78.8 


2.2 


—6.1 


51897.802 


CF /Red 


0.837 


—46.9 


1.3 


1.7 


100.7 


2.2 


-4.4 


51897.910 


CF /Red 


0.840 


—51.4 


1.3 


-0.8 


96.0 


2.2 


— 11.6 


51898.811 


CF /Red 


0.868 


—66.2 


1.3 


0.3 


128.1 


2.2 


-0.8 


51898.922 


CF /Red 


0.872 


—68.7 


1.3 


—0.2 


133.2 


2.2 


1.6 


51899.809 


CF /Red 


0.900 


—83.5 


1.3 


—0.1 


153.0 


2.2 


1.7 


51899.915 


CF /Red 


0.903 


—84.4 


1.3 


0.5 


154.2 


2.2 


0.7 


51900.802 


CF/Red 


0.931 


-95.5 


1.3 


0.6 


171.8 


2.2 


3.5 


51900.908 


CF/Red 


0.934 


-96.3 


1.3 


0.7 


172.6 


2.2 


3.1 


51901.788 


CF/Red 


0.962 


-99.8 


1.3 


0.3 


173.7 


2.2 


0.2 


51901.917 


CF/Red 


0.966 


-100.5 


1.3 


-0.7 


177.2 


2.2 


4.1 


53291.928 


CF/Red 


0.651 


30.8 


2.6 


0.2 








53292.984 


CF/Red 


0.684 


20.6 


0.4 


0.6 








53658.997 a 


2.1m/Blue 


0.187 


52.6 


6.7 


-0.8 








53659.000 a 


2.1m/Blue 


0.187 


51.5 


1.4 


-2.0 








53663.010 a 


2.1m/Blue 


0.313 


71.6 


0.4 


-5.0 








53663.989 a 


2.1m/Blue 


0.344 


89.0 


4.2 


11.7 








53684.903 


CF/Bluc 


0.001 


-85.7 


1.3 


0.8 


155.4 


4.2 


0.1 


53686.902 


CF/Bluc 


0.064 


-30.1 


0.7 


0.9 


76.3 


9.8 


-4.9 


53688.846 


CF/Bluc 


0.125 


21.9 


0.9 


0.9 


13.1 


5.2 


1.3 


54020.972 


CF/Red 


0.564 


51.2 


1.8 


-1.1 


-25.8 


2.2 


4.1 


54024.926 


CF/Red 


0.688 


17.4 


2.8 


-1.0 


13.9 


3.4 


-1.1 


54030.021 


CF/Bluc 


0.848 


-56.0 


2.0 


-0.8 


103.2 


2.4 


-9.3 


54031.999 


CF/Blue 


0.911 


-90.1 


2.9 


-1.7 


144.4 


3.5 


-12.0 
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Table 1 — Continued 



HJD 

(-2,400,000) 


Telescope/ Orbital V\ 
Band Phase (km s _1 ) 


o-i (O - C)i 
(km s _1 ) (km s _1 ) 


V 2 a 2 (O - C) 2 
(km s~ 1 ) (km s~ 1 ) (km s _1 ) 





a zero weight 



Table 2. Orbital Elements for HD 37366 



Element Primary Secondary Joint Solution 



P (days) 31.8187 ±0.0004 31.822 ±0.002 31.8188 a 

Ti (HJD-2,400,000) 53653.013 ± 0.04 ■ • • 53653.02 a 

T 2 (HJD-2,400,000) ••■ 53653.15 ± 0.19 53653.02 a 

ei 0.329 ±0.003 ■•• 0.330 a 

e 2 ••■ 0.35 ±0.012 0.330 a 

tdi (dcg) 211.4 ±0.6 ■•• 211. 6 a 

uj 2 (dcg) ••■ 212 ±2 211. 6 a 

ififkms- 1 ) 88.6 ±0.3 ■•• 88.7 ± 0.2 

^(kms -1 ) ••■ 118.4 ±1.6 117.4 ± 1.2 

^(kms- 1 ) 13.3 ±0.2 ■•• 13.3 ± 0.2 

72 (km s" 1 ) ■•■ 20.6 ±1.3 21.6 ±0.9 

Mi sin 3 i (Mq) 13.8 ± 0.3 ■•• 13.9 ± 0.3 

M 2 sin 3 i (M @ ) ■•■ 10.5 ±0.1 10.42 ± 0.08 

aisini(R ) 52.6 ± 0.2 ■•• 52.62 ±0.13 

a 2 sinj(_R ) ■■■ 69.7 ± 1.0 69.7 ±0.7 

oi (kms^ 1 ) 1.1 ■■• 1.0 

cr 2 (km s- 1 ) ■ ■ ■ 5.3 5.4 



Fixed. 



Table 3. HD 54662 Primary Radial Velocity Measurements 



riJD 

( o Ann nnn^ 
{ — z,4UU,UUU ) 


lelescopc/ 
Band 


Orbital 

Phase 


V r 

(Km s ) 


<r(line— line) 
(km s- 1 ) 


U — O 
(km s~ 1 


46426.830 


CFHT/Ycllow 


0.191 


62.2 




2.0 


-0.3 


46426.904 


CFHT/Ycllow 


0.191 


61.3 




1.6 


— 1.1 


46426.980 


CFHT/Yellow 


0.191 


61.8 




1.6 


—0.6 


46427.802 


CFHT/Ycllow 


0.192 


61.3 




1.7 


— 1.2 


46427.869 


CFHT/Ycllow 


0.192 


61.6 




2.2 


—0.9 


46427.925 


CFHT/Ycllow 


0.193 


61.1 




2.0 


-1.4 


46427.973 


CFHT/Yellow 


0.193 


61.5 




1.6 


— 1.0 


46428.036 


CFHT/Yellow 


0.193 


61.6 




1.4 


—0.9 


46428.132 


CFHT/Yellow 


0.193 


62.0 




2.1 


-0.5 


46428.805 


CFHT/Ycllow 


0.194 


62.2 




2.0 


—0.3 


46429.021 


CFHT/Yellow 


0.194 


61.4 




1.9 


— 1.2 


46429.814 


CFHT/Ycllow 


0.196 


61.8 




1.7 


-0.8 


46429.883 


CFHT/Ycllow 


0.196 


61.5 




2.1 


— 1.1 


46432.853 


CFHT/Yellow 


0.201 


61.4 




2.1 


-1.3 


46432.897 


CFHT/Yellow 


0.201 


60.9 




1.4 


— 1.8 


46432.999 


CFHT/Yellow 


0.202 


62.2 




2.0 


—0.5 


46433.093 


CFHT/Yellow 


0.202 


61.2 




1.5 


— 1.6 


51817.967 


CF /Red 


0.855 


33.2 




0.6 


0.0 


51818.962 


CF/Red 


0.857 


33.5 




2.5 


0.4 


51819.962 


CF/Red 


0.859 


33.7 




1.1 


0.7 


51820.990 


CF /Red 


0.860 


33.0 




2.1 


0.1 


51821.968 


CF/Red 


0.862 


33.0 




0.9 


0.2 


51822.941 


CF /Red 


0.864 


34.6 




3.7 


1.9 


51823.957 


CF /Red 


0.866 


34.5 




3.4 


1.8 


51824.903 


CF /Red 


0.867 


35.2 




1.2 


2.6 


51889.990 


CF /Red 


0.984 


37.5 




2.6 


0.7 


51890.923 


CF/Red 


0.986 


37.1 






0.0 


51892.899 


CF /Red 


0.989 


35.0 




2.0 


—2.6 


51893.926 


CF/Red 


0.991 


37.5 




1.9 


-0.3 


51894.882 


CF/Red 


0.993 


37.5 




1.5 


-0.6 


51894.956 


CF/Red 


0.993 


39.0 




2.3 


0.9 


51895.934 


CF/Red 


0.995 


37.6 




1.3 


-0.8 


51896.033 


CF/Red 


0.995 


39.0 




1.6 


0.6 


51896.881 


CF/Red 


0.996 


39.4 




0.6 


0.8 


51896.952 


CF/Red 


0.997 


37.0 




1.6 


-1.7 


51897.879 


CF/Red 


0.998 


37.9 




1.7 


-1.0 


51897.943 


CF/Red 


0.998 


37.1 




1.6 


-1.8 


51898.891 


CF/Red 


0.000 


39.1 




2.5 


-0.1 


51898.953 


CF/Red 


0.000 


38.3 




2.0 


-0.9 


51899.885 


CF/Red 


0.002 


37.4 




0.8 


-2.1 


51899.947 


CF/Red 


0.002 


38.3 




0.5 


-1.2 


51900.878 


CF/Red 


0.004 


39.0 




2.9 


-0.8 


51900.940 


CF/Red 


0.004 


38.7 




3.1 


-1.1 


51901.885 


CF/Red 


0.005 


39.5 




2.9 


-0.6 


51901.949 


CF/Red 


0.006 


39.1 




2.8 


-1.0 
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Table 3 — Continued 



HJD 


Telescope/ 


Orbital 


V r 


<r(line— line) 


O-C 


(-2,400,000) 


Band 


Phase 


(km s- 1 ) 


(km s- 1 ) 


(km s- 1 ) 


54020.025 


CF/Red 


0.802 


33.9 


7.5 


-2.5 


54024.964 


CF/Red 


0.811 


31.8 


5.6 


-4.0 


54027.025 


CF/Bluc 


0.815 


35.8 


5.0 


0.2 


54028.964 


CF/Bluc 


0.819 


34.4 


6.4 


-0.9 


54030.961 


CF/Bluc 


0.822 


37.4 


4.4 


2.3 


54032.012 


CF/Bluc 


0.824 


35.2 


4.3 


0.2 



Table 4. Preliminary Orbital Elements for HD 54662 



Element 


Value 


P (days) 


557.8 ±0.3 


T (HJD-2,400,000) 


22333 ± 5 


e 


0.28 ±0.04 


w (cleg) 


238 ±5 


K & (km s" 1 ) 


15.9 ±0.5 


7 (km s~ x ) 


49.9 ±0.6 


f{mf (M Q ) 


0.20 ±0.02 


a 1 sini a (i? ) 


168 ±6 


rms (km s _1 ) 


3.3 



a Lower limit due to line blending. 



-19- 



Table 5. Stellar Parameters for HD 37366 



Parameter Value 



V\ sini (km s x ) 
V2 sini (km s -1 ) 
T eSA (kK) 
T efr , 2 (kK) 
log #1 (cgs) 
log 92 (cgs) 
F 2 /F 1 
AM V 



30 ±10 
100 ± 10 
33 ±1 
30 ±1 

4.0 ±0.1 
4.5 ±0.2 

0.35 ±0.05 

1.1 ±0.1 



Table 6. SED Parameters 



Parameter 


HD 37366 


HD 54662 


Primary Type . 


09.5 V a 


06.5 V a 


Secondary Type 


B0-1 V 


09 V 


E(B-V)(mag) 


0.39 ±0.01 


0.32 ±0.01 


R v (mag) 


3.59 ±0.01 


2.82 ±0.01 


9 LD (fi&s) b 


48.4 ±3.0 


72.7 ±3.4 


d (kpc) 


1.38 - 1.92 


1.23 - 1.53 


p max (mas) 


0.4-0.5 


3.7-4.7 



IWalbornl (11972 . 11973 ) 
b Primary 
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ORBITAL PHASE 



Fig. 1. — Calculated radial velocity curves (solid lines) for HD 37366. The primary and 
secondary star's measured radial velocities are indicated by circles (2000), inverted triangles 
(2004 October), stars (2005 October), triangles (2005 November), diamonds (2006 October), 
and squares (IUE 1987). The filled symbols correspond to the primary and the open symbols 
to the secondary. The uncertainties in individual measurements are generally smaller than 
the size of the symbols. 
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Fig. 2. — Tentative radial velocity curve (solid line) for HD 54662 for a period of 558 d. 
The me asured radial velocities a re indicated b y solid circles ( 2000), solid diamonds (2006) , 
squares ( Stickland fc LlovdboOl), op en circles ( Fullerton 1990 ) , stars ( Garmany et al. 1980 ). 



inverted triangles (IConti et al 



19771 ) , and triangles (jPlaskettl 119241 ) . Expanded horizontal 



bars are plotted to show the radial velocities derived from fitting the composite line profile 
from the average spectra for three observational epochs (§5). The dashed and dot-dashed lines 
are the radial velocity curves fit to these time-averaged points for the primary and secondary 
star, respectively. The uncertainties in individual measurements are generally smaller than 
the size of the symbols. 
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0.0 1 , , , i , , , i , , , i , , , I 

6400 6600 6800 7000 7200 

WAVELENGTH (ANGSTROMS) 



Fig. 3. — Tomographic reconstruction of the spectra of HD 37366 based upon the 30 red 
spectra obtained in 2000. This plot shows the primary (top) and the secondary spectrum 
(bottom), as well as absorption line identifications (vertical marks). The atmospheric lines 
in the region of 6850— 7000 A are replaced with the continuum. 



-23- 



2.0| 




0.0 1 I 

4200 4300 4400 4500 4600 

WAVELENGTH (ANGSTROMS) 

Fig. 4. — Tomographic reconstruction of the spectra of HD 37366 based upon five blue 
spectra from runs in 2005 and 2006. This plot shows the primary (top) and the secondary 
spectrum (bottom) as well as a number of absorption line identifications (vertical marks). 




Fig. 5. — Mean spectrum of HD 54662 from our observations. Line identifications are marked 
by vertical lines. 
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Fig. 6.— The mean line profiles of HD 54662 for He I A5876 (left, lFullertonlll990h and He I 
A7065 (right; this work). The observations are plotted as thick lines, the combined model fit 
as thin lines, and the individual component profiles are represented as dotted/ dashed lines 
for the primary and secondary, respectively. 
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1000 10000 

WAVELENGTH (ANGSTROMS) 



Fig. 7. — Best fit SED for HD 37366. The solid line indicates the combined model flux for 
the binary, and the diamonds represent the photometric observations (described in the text). 
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1000 10000 

WAVELENGTH (ANGSTROMS) 



Fig. 8. — Best fit SED for HD 54662. The solid line indicates the combined model flux for 
the binary, and the diamonds represent the photometric observations (described in the text). 



